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We study the non-equilibrium transport properties of a one-dimensional array of dissipative quan- 
tum dots. Using the Keldysh formalism, we show that the dots' dissipative nature leads to a spatial 
variation of the chemical potential, which in disordered arrays, breaks the invariance of the current, 
/, under bias reversal. Moreover, the array's nanoscopic size results in an algebraic low-temperature 
dependence of /. Finally, we show that a local Coulomb interaction splits the dots' electronic levels, 
resulting in a Coulomb blockade, which is softened with increasing dissipation and array size. 

PACS numbers: 73.63.-b, 73.63.Kv, 73.22.-f, 73.23.-b 



Arrays of metallic [T, ^2^ or semiconducting [s*, ^ quan- 
tum dots (QDs) have attracted significant interest due 
to the unprecedented experimental control in assembling 
such arrays, and the resulting ability to custom-design 
their electronic structure and transport properties [5, 6j. 
The continued miniaturization of such arrays, important 
for many applications in quantum computation f?\ and 
nanoelectronics [5] , raises the important question of how 
the non-equilibrium transport properties in nanoscopic 
arrays differ from those in mesoscopic ones and how 
they evolve across different length scales. Nanoscale ar- 
rays also provide a unique opportunity to study how 
the local interplay between dissipation [8, 9], disorder, 
and Coulomb interaction affects their global transport 
properties. Investigating the evolution of these prop- 
erties across different length scales might also hold the 
key to identifying the (yet unresolved) origin of the low- 
temperature conductivity observed in mesoscopic QD ar- 
rays [H, ^ which is ascribed either to some variable-range 
hopping (VRH) [2, 3] or cotunneling [H, Eq| mechanism. 

In this Letter, we address the above questions by 
studying the non-equilibrium transport properties of one- 
dimensional (ID), nanoscale arrays of dissipative QDs, 
as schematically shown in Fig. [T]^a). Using the Keldysh 
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FIG. 1: (a) One-dimensional array of quantum dots coupled 
to two leads, (b) First and (c) second order fermionic self- 
energy corrections due to a local Coulomb interaction. 



Green's function formalism [lll,ll2|, we demonstrate that 
the interplay between dissipation, Coulomb interaction, 
and disorder gives rise to a series of novel non-equilibrium 
quantum effects. In particular, the dots' dissipative na- 
ture leads to a spatial variation of the chemical potential, 
/i, reflecting the array's resistance, which changes qualita- 
tively with increasing disorder, and breaks the invariance 
of the current under bias reversal. Moreover, the array's 
nanoscopic size yields an algebraic low-temperature de- 
pendence of the current in disordered arrays, in contrast 
to the exponential scaling observed in mesoscopic arrays 
[il,|3|. Finally, we demonstrate that a local Coulomb in- 
teraction results in a splitting of a dot's energy level, 
which continuously varies with the level's occupation. 
This splitting gives rise to a Coulomb blockade in the 
charge transport, which is softened with increasing dis- 
sipation and/or size of the array. These results provide 
important new insight into the complex non-equilibrium 
transport properties of nanoscale systems. 

We consider a one-dimensional array of N quantum 
dots connected to a left and right lead [Fig. [Il^a)] whose 
Hamiltonian is given hy H = He + Hph + Hieads where 



~ S>^i+1'^ ~ ^L,adu,a + H.C. (1) 

j,cr a,v=R,L 

Here ^ and dj,^ create an electron with spin a on dot 
j and lead v {v = R^L), respectively. We assume for 
simplicity that on each dot, there is only a single elec- 
tronic state with energy £{j) that is relevant for charge 
transport, as is the case in semiconducting QDs [9]. U 
represents the strength of the local (intra-dot) Coulomb 
interaction, and t is the hopping matrix element between 
neighboring dots, ti, describes the coupling of the array 
to the leads with n/,(^) = 1{N). Finally, Hph and Hieads 
represent the local electron-phonon interaction on each 
dot and the leads' electronic structure, respectively. 

In order to study the array's non-equilibrium transport 
properties, we consider a (symmetric) bias, V = ijll — IJ^r^ 



2 



across the array arising from different chemical potentials 
in the left and right leads with (Ir = —fiL- The resulting 
current between dots j and (j + 1) is given by [13j 
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Here, F is the full Keldysh Green's functions matrix of 
the entire array, which accounts for (electron-phonon and 
Coulomb) interactions, electronic hopping between dots 
and the coupling to the leads. Since scanning tunneling 
spectroscopy (STS) on QD arrays [9] reported that intra- 
dot interactions are larger than the inter-dot electronic 
hopping [9I, we include the former first in the perturba- 
tive calculation of F, in contrast to the limit considered 
previously [HI- Moreover, STS experiments [H, [9| 
also reported significant energy broadening (i.e., dissi- 
pation) of the dots' electronic states (likely arising from 
electron-phonon coupling) which is only weakly tempera- 
ture and energy dependent, and remains substantial even 
for T ^ [9]. We account for the dissipation by intro- 
ducing a (phenomenological) lifetime, F"^, of the dots' 
electronic states (we take £ to include a possible shift due 



to the electron-phonon interaction) [lj,ll5|- The Green's 
function of a single (isolated) dot is then (the effects of 
a Coulomb interaction are discussed below) 



^{uo) = uo - {e - 11) ^ iV . 



(3) 



Including next the electronic hopping between the dots 
and the coupling to the leads, we obtain 



F = (1 - gri) f (1 - iga 



(4) 



Here, gr,a and / are the (diagonal) retarded, advanced 
and Keldysh Greens function matrices describing the de- 
coupled (t = ty = 0) dots and leads, respectively, with 



f{uj) = 2i[l- 2fiF{uj)] Im [gr{uj)] 



(5) 



fip is a diagonal matrix containing the dots' and leads' 
Fermi-distribution functions, and t is the symmetric hop- 
ping matrix. While the results shown below were ob- 
tained in arrays with = 20 dots, we found qualitatively 
similar behavior up to the largest arrays {N = 150) we 
studied. For concreteness, we take t = Ir^l = O.SEq and 
F = 0.05£^o, unless otherwise specified. 

In order to understand an array's transport properties, 
it is instructive to first consider its (equilibrium) elec- 
tronic structure. To this end, we present in Fig. [2fa) the 
local density of states (LDOS), Nj{uj) = —lmG^j{uj)/7r 

where Gr = (l — 9ri) ^ 9r of a non-disordered array with 
N = 20 dots. The array's nanoscopic ID structure gives 
rise to a significant curvature of the band (see discussion 
below), with the LDOS exhibiting the precursors of the 
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FIG. 2: QD array with N = 20. (a) Nj(uj) in equilibrium with 
£(j) = ^(j) = 0. (b) /i(j), s(j) for V = 0.4^0/e and different 
F. L and R refer to the left and right leads, respectively, 
(c) Resistance = 1/// as a function of N . Inset: Rd as a 
function of F. (d) /l/-curve for a non-disordered array and 
two different disorder realizations with s — 0.2Eo. 



1/ V divergence characteristic of an array with N = 00. 

We next discuss the array's non-equilibrium transport 
properties. For a given bias V, charge conservation in 
the steady-state requires that all currents be the 

same. For dissipative QDs with F 7^ 0, this requirement 
can only be satisfied by allowing spatial variations in the 
local chemical potential [16], /i(j), of each dot [Fig.lSfb)], 
which enters in the calculation of through f{uj) in 

Eqs.dl]) and (|5]). Hence, fi should be interpreted as the 
electro-chemical potential and its difference between two 
adjacent dots as the change in the electrostatic potential 
due to the array's resistance. As expected, the spatial 
variation of /i decreases with decreasing F and vanishes 
for a non-dissipative array with F ^ [12| . A plot of the 
resistance, = F//, as a function of N [Fig.lSfc)] reveals 
a linear relationship R = 2Ri-\-R(iN ^ with Rj ^ 0.50h/e^ 
and Rd ~ 0.084fi/e^ being the interface resistance and 
resistance per dot, respectively, demonstrating that the 
array behaves as a series of resistors. Since Rj Rd, 
the change in /i at the interface is significantly larger 
than that between dots, in agreement with Fig. EJb). 
Moreover, Rd ^ T (see inset) demonstrates that it is 
the dissipative nature of the QDs that is responsible for 
the array's resistance. Finally, a plot of the IV-cuive in 
Fig.[2fd) shows ohmic behavior, / ~ at low bias while 
/ becomes basically independent of V once V exceeds the 
electronic band width, W ^ At = 1.2Eo^ of the array. 

Disorder in an array can emerge from non-uniformity 
in the dots' spacing or size, leading to variations in t and 
6:(j), respectively. Below, we consider the latter case, 
which has a more profound effect on the array's trans- 
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FIG. 3: s{j), for (a) a non- disordered array and disorder 
case /, and (b) for a strongly disordered array with t = F = 
0.01^0, and s — 0.5^0- (c) s{j)^ 0) ^ disordered array 
under bias reversal, (d) /(T) for an array with = 150 at 
intermediate disorder with s — t — 0.3^o and F = 0.05^o- 
Inset: /(T) for the strongly disordered array of (b). 



port properties, and take £{j) to be Gaussian (disorder-) 
distributed with average value {e{j)) = and standard 
deviation s. The IV -curve for two realizations of inter- 
mediate disorder (t ~ 5), case I and II with s = 0.2Eo^ are 
plotted in Fig.[2jd) [s{j) for case I is shown in Fig.[3l^a)]. 
While in both cases, the current is suppressed by disor- 
der, the extent of the suppression depends on the specific 
disorder realization. The microscopic origin of this sup- 
pression lies in the diminished hybridization of electronic 
states between neighboring dots due to the differences 
in their energies, s{j). In Fig. [3l^a) we present e{j) and 
fid{j) for disorder case I together with for the non- 
disordered array. Surprisingly, fid{j) does not differ sig- 
nificantly from and in particular, does not track the 
variations in s{j). The reason for this robustness of jiid 
is that in the intermediate disorder case, electrons can 
tunnel even through those dots where \e\ is large. In 
contrast, in the strong disorder limit, t, F <C 5, shown in 
Fig.El^d), fid{j) exhibits step-like changes between those 
neighboring dots where the variation in e{j) is the largest. 
This is expected since large spatial variations of e{j) in- 
dicate the regions of largest resistance in the array, with 
concomitant large (potential) changes in fid{j)- 

The interplay between disorder and dissipation gives 
rise to an interesting novel quantum phenomenon in 
which the invariance of the current's magnitude is bro- 
ken under bias reversal. To demonstrate this effect, 
we present in Fig. [Sjc) the chemical potentials, /ij(j), 
in a dissipative, disordered array for two different bias 
V± = :^Eo/e. Due to the combination of disorder and 
dissipation, /ij" and /i^ are not related by a spatial sym- 



metry, implying that the corresponding currents are in 
general different: in Fig. [3fc), the current for positive 
bias is /+ = 0.0089e£^o/^, while for V- one finds 
/_ = —0.0261eEo/h. We note that the invariance is re- 
stored either in a non-disordered dissipative array, or in 
a non-dissipative (F = 0+) disordered array. In the latter 
case, the current is independent of /i and /ij(j) = 
again satisfies the spatial symmetry. 

The temperature dependence of / is generally deter- 
mined by both disorder and interaction effects, with the 
latter being represented by F. However, when disorder 
becomes stronger than dissipation (i.e., s > F), we find 
that the current is predominantly determined by the for- 
mer, and becomes essentially independent of F. Thus, in 
calculating /(T), we can neglect the temperature depen- 
dence of interaction effects, i.e., r(T), as was previously 
assumed in Refs.[l3, Ell. In the limit of intermediate 
disorder, F <C t ~ 5 , and for ksT <C s, I{T) possesses 
an algebraic temperature dependence, /(T) = /q + BT^^ 
as shown in Fig. [Sjd) for an array with N = 150 dots. 
This algebraic dependence is the lowest order (in T) con- 
tribution arising from the frequency dependence of the 
Green's functions entering Eq.(j4]), which in turn arises 
from the finite ID structure of the array [this frequency 
dependence is also present in the LDOS of Fig. [2fa)]. In 
the strong disorder limit, F, t ^ s (see inset of Fig. [31(d)], 
the algebraic temperature dependence of /(T) at small 
T is followed by a more rapid increase at larger tem- 
peratures, which can be described by an activated be- 
havior over some limited temperature range. We note 
that the absence of an exponential VRH temperature 
scaling 17, ill, is one of the characteristic hallmarks of 
nanoscale systems. Not only would the algebraic con- 
tribution to I{T) dominate any exponential scaling for 
T ^ 0, but VRH scaling also assumes that a disordered 
system is characterized by a single parameter, the local- 
ization length, and thus neglects spatial fluctuations in 
the disorder. These fluctuations are, however, highly rel- 
evant in nanoscale systems due to the lack of disorder self- 
averaging, resulting in different /V-curves for different 
disorder realizations with the same s [Fig. [21(d)]. Work is 
currently under way to identify the characteristic size of 
the QD array at which disorder self- averaging will lead 
to a crossover from algebraic to VRH scaling of I{T). 

We next study the effects of a local Coulomb interac- 
tion by including the first [Fig. [11(b)] and second order 
[Fig. [11(c)] fermionic self-energy corrections in the local 
Green's functions entering Eq.(j4]). Based on the experi- 
mental STS results of Ref. [9] , we take the Coulomb inter- 
action to be weaker than the electron-phonon interaction, 
and thus take the internal fermionic lines of the diagrams 
in Figs.[Tl(b) and (c) to be given by the dissipative Green's 
function of Eq.(|3]). The first order diagram of Fig. [Tl(b) 
leads to an energy shift ea{j) = £{j) + Unj^a- with rij^a- 
being the electron occupation and cr, a being opposite 
spin quantum numbers. In contrast, the second order di- 
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FIG. 4: (a) LDOS of a single dot with F = 0.05^o for = 
0, 0.4, 0.5. (b) IV-cuTve for an array witli N — 1 and different 
r. (c) Normalized IV-cuive for arrays with N — 1 and N — 5 
dots, (d) A(j) = — s{j) as a function of V for N = 5. 



agram, for which an analytical form can be derived, leads 
to an energy splitting which depends on rij^a' for a com- 
pletely empty or filled state, the splitting vanishes, while 
it is the largest for a half- filled state. To demonstrate 
the evolution of the dot's electronic structure with n^^cr, 
we present in Fig. Hl^a) the LDOS of a single dot [in- 
cluding the self-energy corrections of Figs. [Hb),(c)] for 
= 0.0,0.4,0.5, and U = 2Eq (the LDOS for occu- 
pation 1 — rij^a- is obtained via uo ^ —uj). For Uj^o- = 
(^j,cr = 1)7 the dot possess a single, doubly degenerate en- 
ergy level at cj = —Eq {u = -\-Eo). As rij^a- increases, this 
level shifts to larger energies, and a splitting emerges. For 
^j,cr =0.5, the state's spectral weight is predominantly 
located at = ±£^0, i-e., separated hy U = 2Eo^ while a 
small contribution of spectral weight is found at = 0, 
which is a direct consequence of F 7^ 0, and hence van- 
ishes for F ^ 0. This evolution of the LDOS with rij^a- is 
identical to that of Ref. [19] for F = 0+, but qualitatively 
differs from it for F 7^ 0. The state's splitting, which 
is analogous to the Coulomb interaction driven suppres- 
sion of the LDOS near the Fermi energy discussed by 
Efros and Shklovskii [l8|, is the origin of the Coulomb 
blockade observed in IV curves. For a single QD with 
/i = and rij^a- = 0.5 [Fig. ID^b)], the dot's energy level 
is maximally split, and / is strongly suppressed until V 
exceeds [/, thus manifesting the Coulomb blockade. The 
state's spectral weight at = (due to F / 0) leads to 
a non-zero / even for a V < U. The Coulomb blockade 
is softened both with increasing F [Fig. ID^b)] as well as 
with increasing number of dots N in the array, as shown 
in Fig. m^c) for = 5. In addition, in the latter case, 
new features appear for V < U which arise from changes 
in A(jf) = /i(jf) — £(j), where A(j) = corresponds to 



^j,cr = 0.5 and thus a maximum splitting of the energy 
level. As A(jf) increases, this splitting is reduced, and / 
increases (and vice versa) as follows from a comparison 



of Figs. lU^c) and (d). Therefore, the local maximum of 
/ at Vi ~ 0.7 Eo/e coincides with that of A(jf), while its 
sharp increase at V2 ^ 1.2Eo reflects that of A(j). For 
V > I decreases slightly due to an increasing separa- 
tion between the dots' energy levels, but then increases 
again as V approaches the Coulomb gap, U. 

In summary, we studied the non-equilibrium transport 
properties of a ID array of dissipative QDs. We showed 
that the QDs' dissipative nature leads to a spatial varia- 
tion of /i(j), reflecting the resistance of the array. In dis- 
ordered arrays, this variation breaks the invar iance of the 
current under bias reversal while the array's nanoscopic 
size results in an algebraic low-temperature dependence 
of I{T). Finally, we showed that a local Coulomb interac- 
tion gives rise to a splitting of the dots' energy levels, and 
the emergence of a Coulomb blockade, which is softened 
with increasing dissipation and size of the array. 
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